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1 Summary 

On Ist Jan 2014, the rate of production in the area of Loppersum was reduced. Here we seek to understand if 
reduction in production has changed the characteristics of inter-event times for seismic events in the vicinity 
of Loppersum. This reports shows that if only events of magnitude >= 1.5 are considered, there is little 
or no evidence for change in production having influenced inter-event time. If all events (i.e. events of all 
magnitudes) are considered (as reported in the KNMI catalogue), there is considerable evidence that the 
so-called ’catalogue inter-event time’ after production change is different to the earthquake inter-event time 
immediately prior to the production change. Specifically, we find that events of small magnitude were more 
frequent before production change than after. This observation should be viewed with considerable caution, 
since (1) recording of events of magnitude < 1.5 is uncertain both in spatial and time domain (according to 
the information provided by KNMI), and (2) the period of observation since production change is short. 

It is further recommended that a further study of inter-event time with magnitudes >= 1.0 should be 
conducted to understand the influence of change in production more fully, for carefully-specified spatial 
subdomains. These spatial subdomains should include locations within the Loppersum boundary where 
trusted geophones are located, partially reducing uncertainty in our analysis for events of magnitude <1.5 
(since we do not know the exact location of the geophones). 


2 Methodology 

2.1 Inter-event time 


Let there be n number of events after the production was reduced on Ist Jan 2014 up to the current date 
(18th Jan 2014, the last date in the event catalog) in the area defined by the production clusters covering 

Loppersum. Let ti, t 2 , . be the inter-event time from a sampling density f(ti,t 2 , - JnA) where A is the 

parameter and one assigns A a prior p{X). Bayesian posterior density of A is given by the following equation 


p{X\ti,t2...tn) 


f{tl,t2,-tn\X)p{X) 

f{ti,t2,..tn\X)p{X)dX 


The above equation can be simplified to proportionality so that: 


( 1 ) 


p{X\ti,t2...tn) OC f{ti,t2, ..tn\X)p{X) (2) 

If activity rate follows a Poisson process (ref Bourne and Oates Activity rate model), ti, t 2 ■■ tn are from 
Exponential distribution. 


f{tl,t2,-tn\X) = /(t,|A) OC 


(3) 


I 
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This can be written as: 

--inlA) (X A“”e^ (4) 

where s is the summation of the n inter-event time. 

Substituting 4 in 2 gives the posterior density of A as: 


p{X\ti,t 2 .:t„) (X X ”e ^"p(A) (5) 

The prior density function for A is unknown and one can assume a non informative prior up to proportionality 
as: 

p(A) oc j (6) 

which when substituted in 5 leads to the posterior density 

p(A|ti,< 2 --d„) oc — (7) 

or 

p(A|ti,t 2 ---in) oc (8) 

It is often beneficial to deal with log-likelihood for the posterior density in order to deal with small numbers 
which leads to: 


log{p{X\ti,t 2 ...tn)) oc —(n + l)logX 


s 

X 


(9) 


It is also possible to assume that the prior density p{X) is the non-informative Jeffrys’ prior which would lead 
to the posterior density to be slightly different. 


log{p{X\ti,t 2 ...tn)) oc -(n + -)logX 


s 

A 


( 10 ) 


Also, one can assume the prior density to be uniformly distributed which leads to 

g 

log{p{X\ti,t 2 ...tn)) oc -nlogX- - (11) 

It should be clear from equations 9, 10 and 11 that the assumption about the prior can influence the pos¬ 
terior density to some extent. Also, for an exponential variate, an informative prior can come from gamma 
distribution with two known non-zero shape parameters a and /3. While we can estimate them for the period 
prior to production change, it will be hard to estimate them for the period after the production was changed 
due to the limited sample size. For the purpose of the this report, we will therefore use equation 9 while 
posteriors based on other non informative or informative prior can be computed at a later stage if required. 
Also, for the saké of simplicity, we henceforth refer to the posterior density p{X\ti,t 2 ---tn) simply as p{X). 

In order to estimate the posterior density, we follow Markov Chain Monte Carlo (MCMC) approach. To 
do so, we initiate a sequence of A with an initial value A*^ which could be a random draw from a uniform 
distribution. In the next step, we simulate a candidate value A* as the 1®* proposal in the sequence from a 
normal proposal density (/(A*|A*^) and compute the difference 


R = log{p{X*)) - log{p{X°)) 


(12) 
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If i? > log{u) where u is a simulated random uniform random variate [0, 1], then the proposal A* is accepted 
in the sequence and is A* otherwise A^ is A°. We now repeat the steps i.e. generate A*, from normal 
proposal density g(A*|A^) and compute the difference : 

R = log{p{X*)) - log{p{X^)) (13) 

A^ is A* if i? > log{u) otherwise A^ is A^. In this way, A^, .... A^ where k is some large number are simulated 
in the sequence. The end of the chain is stable and can be used to compute the kernei density of A. 

The above discussion for estimating the posterior density for A after the production change occurred can 
also be used to estimate the posterior density of A before the production change between the period 2003 to 
Dec 2013. In this case, if N is the total number of events for the entire period (2003 to 2013) but in the 
same area that defines n, N »> n simply because of the time period prior to the production change has 
allowed large number of events to occur over time. So, to estimate A prior to the change in production in 
such a way so that it can be compared with A after the change in production, we simply divide the N events 
into k contiguous blocks, each of size n and use equation 9 for each block. There will be k A’s and Xk -2 to 
Afc should be informative for comparison with A after the production was changed on Ist Jan 2014. This is 
because A for blocks k — 2 to k which refer to the time period before the production change are closest in 
time to A for the block after the production was reduced and the influence of other covariates on the seismic 
activity rate can be ruled out. If the density functions of A before and after the production change are over¬ 
lapping, we can conclude that the weight of evidence suggests that production change has resulted in a small 
change in the activity rate up to the current date. On the other hand if the density functions do not overlap, 
we can infer that the weight of evidence for a change in activity rate due to change in production rate is large. 


2.2 Defining the Influence Area 

We now need a systematic way to define the area whose seismic activity can be influenced by production 
change in Loppersum as this defines n , N and s used to estimate A. In all there are 22 production clusters 
in the Groningen field. Five of these production clusters namely Leermens, Overschild, Ten Post, t Zandt 
and De Pauwen belong to Loppersum. They are north west of Groningen field and their production was cut 
on Ist January 2014 as shown in figure 1. Raw production is on a daily basis and smooth production is a 
30 day moving average with zero phase lag. Notice that the production post Jan 14 does not remain at zero 
at all times and it seems some volume is produced during high demand, winter season. Also, in pre-2014 
period, the production for each cluster is highly cyclical with peaks near December/January and troughs in 
a summer month. 

There are three possible ways to understand the seismic area influenced by these five clusters qualitatively. 
This is shown in figure 2. In the first scenario, the area of seismic influence is assumed to not only cover the 
interior area drawn by the polygon that defines the five clusters but also extends in the north west direction. 
This is a fair assumption as there are no production clusters north west of Ten Post and t Zandt. It is 
therefore safe to say that only production from these two clusters can cause seismic events in north west 
direction if we assume a causal relationship between production and events. In the second scenario, the area 
of seismic influence is assumed to be moderate and extends only few kilometers to the north west. In the 
third scenario, the area within the polygon and a small area outside the polygon is assumed to be influenced. 
Notice that in all the three scenarios, the area influenced is not overlapping with other production clusters. 
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Figure 1: Raw and smoothed production by clusters belonging to Loppersum and their sum (bottom right) 
for the period 2003 to Jan 2015. 





















































































DRAFT NOTE 





O All Other clusters 

# Leermens 

# Overschild 
O TenPost 

# tZandt 

# De Pauwen 

— Loppersum boundary 
Area of Seismic influence 


Figure 2: Clusters belonging to Loppersum coloured circles. Production changes to these clusters is assumed 
to influence seismic activity in the area covered by cyan rectangles of different sizes. Three scenarios are 
considered. Top left: Large influence area Top Right : Moderate influence area. Bottom left: Narrow 
influence area. Notice influence areas are not overlapping with other clusters. 

2.3 Treating for aftershocks 

A seismic event can also trigger another seismic event. This is often called as an aftershock. It is important 
to understand that an aftershock is not a true event but just a consequence of another event. An aftershock 
event has to be properly identified and removed from our analysis as it can give a biased view of the the 
underlying inter-event time. There are various methods suggested in the literature to identify an aftershock 
but for the purpose of our work, we have treated any event which occurs within three days of the preceding 
event as an aftershock. Thus, by this methodology, the minimum inter-event time in any block has to be 3 
days. 


3 Events of Magnitude >=1.5 

We now try to understand if the change in production has an influence on seismicity using inter-event time 
(A) as a metric described in the preceding section on methodology under three possible. In other words, does 
inter-event time A before and after production change when the area influenced by the production change is 
assumed to be broad, moderate and narrow? 

3.1 Broad influence area 

When the seismic area influenced by the change in production in Loppersum clusters is assumed to be broad, 
the spatial distribution of seismic events is shown in figure 3. There are (n =) 8 events of magnitude >= 1.5 
after the production was changed in Jan 2014. The events that occurred prior to Dec 2013 can be arranged 
into (k =) 11 blocks of 8 events each such that the last three blocks (blocks 9, 10 and 11) are closest to 2014. 
Figure 4 shows the temporal distribution along with the sum of production for the five clusters while figure 5 
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shows inter-event time, which is also shown in a tabular form in table 1. Posterior A’s, given the information, 
can now be estimated for all the blocks using Markov Chain Monte Carlo (MCMC) technique and their 
smooth kernei density is shown in figure 6. A density for blocks 9, 10 and 11 corresponds to inter-event time 
before the production change was effected while block 12 corresponds to A for the post production change 
period. Figure 7 shows that all density curves of A overlap suggesting that null hypothesis cannot be rejected 
and the available data suggest that change in production has not impacted the activity rate. Figure 8 shows 
the 95 percent confidence interval of A for all the blocks. We can potentially argue that the trend for median 
A is declining with time (activity rate is potentially increasing). However, the spread is large and the median 
for the last block (refers to period after production change) is not any different to the median from all the 
other blocks. Effect of production change, if any, is not evident in this limited time frame. 


.|q 5 Production Clusters, Loppersum boundary, Area of Study, Events Pre 2014 and Post 2014 
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Figure 3: 


Clusters with reduced production. Seismic events before and after production cut. Broad influence 


area. 
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Figure 4: Production(top) and Activity(bottom) with broad influence area. 
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Figure 5: Production(top) and earthquake inter-event time (bottom) before and after production cut with 
broad influence area. 


Table 1: Table of earthquake inter-event time: Broad influence area. Block 12 refers to inter-event time after 
production was changed. Blocks 1 to 11 refer to inter-event time before production was changed. 
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Figure 6: Kernei density for A for each block. Magenta (bottom right curve) is for earthquake events after 
production was changed. The three density curve in Red, Green and Cyan are for three periods prior to 
production change and are in period 2012-2013. 
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Figure 7: Kernei density for last three blocks prior to production change in red, cyan and green. Density 
curve in magenta is for period after the production was changed. AU density curves overlap. 
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Figure 8: 95 percent confidence interval A for all blocks. 

3.2 Moderate influence area 

When the seismic area influenced by the change in production is assumed to be moderate, the spatial 
distribution of seismic events is shown in figure 9. There are (n =) 6 events after the production was 
changed. The events that occurred prior to Dec 2013 can be arranged in(fc =) 13 blocks of 6 each in a 
similar way as discussed in the previous chapter. Figure 10 shows the temporal distribution, figure 11 shows 
inter-event time which is also tabulated in table 2. Posterior A’s are in figures 12 and 13 which shows that 
the densities of A before and after production change overlap. This again suggests that the null hypothesis 
cannot be rejected and the available data suggest that change in production has not impacted the activity 
rate. Figure 14 shows the trend in median A and 95 percent confidence interval. Again, it is difScult to claim 
that median A had changed significantly after production change. 
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Figure 9: Clusters with reduced production, Seismic events before and after production cut. Intermediate 
influence area. 
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Figure 10: Production(top) and Activity(bottom) with intermediate influence area. 
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Figure 11: Production(top) and earth quake inter-event time (bottom) before and alter production change 
with narrow influence area. 


Table 2: Table of earthquake inter-event time: Intermediate Influence Area. Block 14 refers to period after 
production was changed. Block 1 to 3 refer to period before production was changed. 
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Figure 12: Kernei density for A for each block. Magenta (bottom right curve) is for earthquake events after 
production was changed. The three density curve in Red, Green and Cyan are for three periods prior to 
production change and are in period 2012-2013. 
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Figure 13: Kernei density for last three blocks prior to production change in red, cyan and green. Density 
curve in magenta is for period after the production was changed. AU density curves overlap. 












DRAFT NOTE 



Figure 14: 95 percent confidence interval A for all blocks. 


4 Narrow influence area 

When the seismic area influenced by the change in production is assumed to be restricted to the area of the 
polygon and its vicinity, there are only 5 seismic events in the period after the production change. Figure 
15 shows spatial distribution of events, figure 16 shows temporal distribution, figure 17 shows distribution 
of inter-event time while ügures 18 and 19 show posterior densities for A. Table 3 shows inter-event time. 
We again see overlap between density curves in cyan (post production change) and pre-production change 
period suggesting the available data cannot be used to suggest that changes to production have influenced 
the seismic activity. 
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Figure 15: Clusters with reduced production, Seismic events before and alter production change. 
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Figure 16: Production(top) and Activity(bottom), narrow influence area. 
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Figure 17: Production(top) and earthquake inter-event time (bottom) before and alter production change, 
narrow influence area. 


Table 3: Table of earthquake inter-event time: Narrow Influence Area. Block 13 refers to period alter 
production was changed. Blocks 1-12 refer to period before production was changed. 
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Figure 18: Kernei density for A for each block. Magenta (bottom right curve) is for events alter production 
was changed. The three density curve in Red, Green and Cyan are for three periods prior to production 
change and are in period 2012-2013. 
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Figure 19: Kernei density for last three blocks prior to production change in red, cyan and green. Density 
curve in magenta is for period after the production was changed. AU density curves overlap. 
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Figure 20: 95 percent confidence interval A for all blocks. 


5 All events 

In the previous section, we used events of magnitude >=1.5 to understand the relation between production 
and seismicity. In this section, we do a similar analysis but include all the events in the catalog that exists 
on the KNMI website. It should be mentioned that though this analysis is statistically robust, the results 
cannot be treated as reliable since events of magnitude <1.5 cannot be recorded accurately. The inter-event 
time should therefore be read as Catalogue inter-event time throughout the text 

5.1 Broad Influence Area 

Figure 21 shows the spatial distribution of all events that occurred prior to production change and also after 
production change in the area which we believe can be influenced by production. The number of events is 
much greater than when only events of magnitude 1.5 and above are considered. After the production was 
changed, 22 events of small and large magnitude occurred in the defined area. The temporal distribution 
and inter-event time is shown in figures 22 and 23 respectively. We can observe that the events get denser 
near 2011 and continue to remain so immediately after the production change was made . The inter-event 
time for block 11 which falls in post production change period is trending upwards as seen in Figure 24 which 
implies that production change has perhaps reduced the number of seismic events of smaller magnitude after 
production change. However, we must wait for a longer period of time in order to confirm this trend. 
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Figure 21: Clusters with reduced production, catalogue events before and alter production change. 
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Figure 22: Production and catalogue events (Yes = 1, No = 0). 
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Figure 23: Production and catalogue Inter event time. 
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Figure 24: 95 percent confidence interval of A. 


5.2 Moderate Influence Area 

Figure 25 shows the spatial distribution of the seismic events that have occurred in the area of influence 
before and after production change was made. Figure 26 shows the temporal distribution of events. We 
again observe that the events occurred more frequently from 2010 to end of 2013, the date of production 
change. After the production change, they were relatively less frequent as seen in figure 27. Figure 28 shows 
that there is a good evidence to believe that the inter-event time has increased after production change but 
again, we must wait for a longer period before confirming this observation. 
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Figure 25: Clusters with reduced production, catalogue events before and after production change. 
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Figure 26: Production and catalogue events (Yes = 1, No = 0). 
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Figure 27: Production and catalogue inter-event time. 
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Figure 28: 95 percent confidence interval of A. 

5.3 Narrow Influence Area 

Figure 29 shows the spatial distribution of all the events in what we believe to be a narrow influence area 
after production change was made. We again see that the temporal distribution of events became dense 
since 2011 (Figure 30) in this area and continued to remain so just after production change was made. Later 
during the year, events seem to get less dense which is also supported by Figure 31 and Figure 32. 
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Figure 29: Clusters with reduced production, catalogue events before and after production change. 










Production mscf/d 


DRAFT NOTE 



2015 


Production in ciusters ciose to Loppersum 


2004 


2006 


2009 


2012 


Events in Area of seismic study 



Figure 30: Production and catalogue events (Yes = 1, No = 0). 
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Figure 31: Production and catalogue inter-event time. 
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Figure 32: 95 percent confidence interval of A. 


6 Remarks and Notes 

The current study shows that the inter-event time of seismic events of magnitude >=1.5 after the production 
change was made has not changed but it appears that if all the events are considered, the seismic activity is 
reduced under all possible scenarios of seismic area influence influenced by production. We should, however, 
caution that though the statistical approach and results are robust, the study is covering a limited time 
period (January 14 to Jan 15) with reduced production. The sample size, in terms of number of events in 
this period is small which can influence the outcome. Also, it is uncertain if events of magnitude < 1.5 
can be recorded accurately and there can be spatial and temporal uncertainty in the catalogue provided by 
KNMI website. The inter-event time when all events are considered should therefore be read as Catalogue 
inter-event time. 

Production for the clusters belonging to Loppersum (and for other clusters as well) is highly cyclical. Each 
year, the peak in production occurs in one of the winter months and trough in summer. Despite the troughs 
in summer, it is not the case that summer months are devoid of any activity. Also, while Dec and Jan months 
have peak production, it should be kept in mind that December is the month of lowest seismic activity as only 
5 seismic events have occurred in this month over 2003 to 2013 period. On the other hand, January shows 
high seismic activity and at-least one event has occurred in this month for all the years from 2003 to 2013. 
This strange observation shows that other covariates such as seasonality in seismicity besides production are 
(also) involved and should be given due importance. 

In the present work, we have assumed that the posterior A is an Exponential variate with non-informative 
prior that is inversely proportional to A. Eor exponential distribution, an informative prior can be based on 
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Gamma distribution with two hyper-parameters. If we have prior belief about the hyper-parameters, the 
analysis can be re-run although we believe there will be a minimum impact on the analysis and conclusions. 

Inter-event time in just one metric by which we can understand the influence of production change on 
seismicity. Other metrics are total seismic moment, frequency of large magnitude events and correlation 
between production and events. These need to be considered to have a proper view on the relationship 
between production changes and seismicity. 

It is recommended that a further study of inter-event time of a spatial subset of events of magnitude >= 1.0 
should be conducted to understand the influence of production. This spatial subset should come from loca- 
tions within the Loppersum boundary where the geophones are present. This will overcome the uncertainty 
in our analysis for events of magnitude < 1.5 as we do not know the exact location of the geophones used to 
record seismic events. 



